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I. INTRODUCTION 

The control of quantum systems and their applica- 
tions have rejuvenated the study of natural and artificial 
atomic systems. One of the most exciting such studies is 
the light-matter interaction, exemplified by cavity-QED 
and the Jaynes-Cummings Hamiltonian pQ. In particu- 
lar, the use of superconducting resonators, with harmonic 
oscillator modes of the electromagnetic field in the mi- 
crowave domain, has gone through a remarkable trans- 
formation in the past decade. 

Early on, harmonic oscillator modes were incorporated 
in superconducting qubit circuits as way to couple them 
together, and thus perform joint operations. The first 
such proposal [5] used a dispersive coupling in which the 
qubit frequency, uj q , is much less than the resonator fre- 
quency, uj r , and this design was refined in [3]. Reso- 
nant coupling was first discussed in [3], and yields an 
architecture that is very similar to the ion-trap quan- 
tum computer pj], m which the resonator is the ana- 
logue of the center-of-mass mode of the linear ion chain. 
The qubit 's quantum state can be transferred to the os- 
cillator, which then interacts with other qubits, and is 
finally transferred back. Resonant coupling was later 
adapted to couple charge qubits [5] and flux qubits [7] 
via an LC-oscillator, and charge qubits via a current- 
biased junction 0[9]. In parallel with these fundamental 
studies came the observation that these systems provide 
strong analogies with cavity quantum electrodynamics 
(QED) [22 HI]- This general approach has now come 
to be called circuit-QED [12], and has led to impressive 
experimental progress [T3H2T] . 

The use of on-chip superconducting resonators to mea- 
sure and couple qubits was then extended to include ideas 
for producing a quantum memory element and bus for 
transferring information, with either microwave |171 122] 
or even nanomechanical [2"3"H2l)] modes. This has culmi- 
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nated in the "von Neumann" model of quantum logic by 
Mariantoni et al. [2BJ, consisting of a qubit-based pro- 
cessing unit, a resonator-based memory, and a resonator 
data bus. However, while the nonlinearity of a qubit or 
other auxiliary is needed to access the individual states 
of the resonator, there remain alternatives to qubit-based 
logic. 

Quantum systems with d > 2 logical levels are tra- 
ditionally called qudits, and unitary gate synthesis [27] 
and error correction |28j can be constructed using qu- 
dits. The Fock states of a superconducting resonator can 
be addressed as a qudit [29], with the potential to out- 
perform qubit-based processing. It is the purpose of this 
paper to study one of the simplest such tasks in quantum 
information processing, namely the generation of entan- 
glement between two systems. While there is essentially 
just one type of entangled state between two qubits, qu- 
dits can be entangled in any of the states 

N 

|*) = ^c n |n)®|iV-n), (1) 

with an entanglement up to \og 2 (N + 1) ebits [30]. We 
shall study the synthesis of the maximally entangled two- 
qudit state (with c„ = 1/y/N + 1) and the "NOON" 
state 

|*noon> = "^(PV)® |0> + |0)®|iV)), (2) 

a state of interest for quantum metrology [31] and re- 
cently demonstrated [32] in a superconducting circuit. 

The NOON state experiment followed the implemen- 
tation [33 of the quantum state synthesis algorithm for a 
single oscillator coupled to a qubit proposed by Law and 
Eberly [34j . In this implementation, a superconducting 
phase qubit was capacitively coupled to a superconduct- 
ing resonator. By performing a sequence of "shift pulses" 
and Rabi pulses, the qubit was put into resonance with 
the resonator, thus transferring quanta between the two 
systems. It was shown in [34 that a properly chosen 
sequence of such operations could generate an arbitrary 
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state of the resonator. In the experiment of Hofheinz et 
al., several quantum states were prepared and analyzed 
by Wigner function tomography 33J . An extension of 
this method, using two three-level systems coupled to 
two resonators [55] , was used to generate the NOON state 
with N = 3 [32] . 

An alternative approach was previously proposed by 
the authors [36 to synthesize an arbitrary entangled state 
of two resonators. Our algorithm was based on previous 
studies of Law-Eberly-like schemes to generate two-mode 
quantum states of the motional degrees of freedom of a 
trapped ion [57J IM| ; suitably modified to use the inter- 
actions available in superconducting circuits. In particu- 
lar, the number-state-dependent transitions in the quasi- 
dispersive regime demonstrated by Johnson et al. [59"] 
played a critical role in allowing selective manipulations 
of the Fock states of the two resonators. 

An important issue missing from previous analyses 
[551 155] is the role of decoherence in these methods. The 
coherence of the coupled system is influenced by deco- 
herence in both the qubit and the resonator, and under- 
standing these effects is an important theoretical ques- 
tion. In this paper we provide a complete analysis of the 
state synthesis algorithm and its application to the entan- 
gled states mentioned above. We will further discuss the 
experimental issues, including decoherence, that arise in 
the synthesis of NOON states, and compare the fidelity of 
the results of the state-synthesis algorithm [36] with the 
NOON state sequence of [35] . The methods and results 
of this study provide a solid base for understanding more 
complex manipulations of superconducting resonators as 
qudits [29] , 

This paper is organized as follows. In Section II the 
state-synthesis algorithm is explored in some detail, with 
special attention given to physical implementation. In 
Section III, we provide a numerical simulation of NOON- 
state synthesis in the strong-coupling regime, illustrating 
some of the subtleties of implementation. In Section IV, 
we present results of NOON-state synthesis in the pres- 
ence of decoherence. Finally, we conclude in Section V, 
while details of decoherence calculations are saved for an 
Appendix. 



II. STATE SYNTHESIS ALGORITHM 

A. Physical Model 

Our algorithm is designed with the following model 
Hamiltonian 

H=H q (t)+Hr+H a +H b (3) 



with 

n q (t) = fiw,(<)o- + o-_ + ifi[n(t)o- + + n*(t)o-_] ) (4) 

W r = huj a a) a + huj^b, (5) 
U a = hg a (a + a + a-a^) , (6) 
H b = hg b (a+b + a-tf) , (7) 

where <r+ = |l)(0|,a^, and b^ are the creation operators 
for the qubit, resonator a, and resonator b, respectively. 
Our goal is to control an arbitrary state of the resonator, 
implementing the transformation 

C/|0,0,0) = £ \Q,n a ,n b ), (8) 

where \q, n a , n b ) = \q) ® \n a ) ® \n b ) is the state in which 
the qubit is in state q — or 1, resonator a is in Fock 
state n a , and resonator b is in Fock state n b . All controls 
are implemented by manipulating the qubit by shifts of 
the qubit frequency uj q (t) or application of resonant Rabi 
pulses of the form fl(t) = n e lUdt . We will shift the 
qubit frequency from the dispersive regime, with ui a < 
Lj q < co b (generalizations will be discussed later), to the 
resonant regimes ui q = ui a or ui q — ui b . In this section we 
will present the theoretical background and motivation 
for the state-synthesis algorithm. An explicit application 
will be presented in the following section. 

In the dispersive regime g a <C | uj q — uj a | , g b -C \uj q — uj b \, 
the states \q, n ai n b ) are approximate eigenstates. The 
Hamiltonian couples each such state to two others (i.e. 
\0,n a ,nb) -t |l,n a -l,n b ), |l,n Q ,n 6 -l) and \l,n a ,n b ) ->■ 
|0, n a + 1, n b ), |0, n ai n b + 1)). Perturbation theory yields 
the following approximation for the energies E qyna ^„ b : 

g 2 

/h « w Q 7i a + u! b n b H n a 

LO a - UJ q 

+ -^n bl (9) 

LU b - UJ q 

g 2 

Ei,n a , nb /ft ~ + u) a n a + to b n b H — — (n a + 1) 

UJ q - LU a 

+ — ^— (n 6 + l). (10) 

LU q - U) b 

The drive frequency needed for the transition 
\0,n a ,n b ) -> |1, n a ,n b ) is thus given by 

= (£"1, — Eo,n a ,n b )/ti 

2 2 

w LOq + (2n a + 1) 9a + (2n b + 1) 

u q — UJ a u q — u b 

(11) 

This Stark shift can be used to selectively address 
individual states (n a ,n b ), conveniently illustrated by a 
Fock-state diagram as in Fig. [T] Here we have set 
5a/( w <? ~ u o) = g 2 / ("ft - u q ) = Aui/2 [54 , allowing us to 
simplify the drive frequencies to 

uj n = ui q + nAuj. (12) 
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For this drive frequency, only those states for which 
n a — n b — n is constant will undergo the transition 
|0, n a ,n b ) — » |1, n a ,nb)- Note that this assumes a two- 
level system for which uj a < ui q < ui b . As discussed be- 
low, this condition is not strictly necessary for the state- 
synthesis algorithm. The key requirement for an efficient 
state synthesis approach is to have an interaction such 
that a given Fock-state (n a ,n b ) can be addressed inde- 
pendently of the neighboring states (n a ± l,n& ± 1) [3"8j . 

We label the Ud — u) n transition RqI . Letting [&) be 
given by 



B. Algorithm 

These three interactions, RqI , 
trated in Fig. [I] By manipulating the qubit frequencies 
and Rabi drives, these three interactions can be effec- 
tively turned on and off sequentially for varying dura- 
tions of time, directing population around the Fock-state 
diagram. By starting from |0,0,0), the target state 



E/|0,0,0>= J2 



t \0,n a ,n b ) 



(20) 



I*) 



c q,n a 

q,n a ,n b 



the transition has the following effect: 

R ot\^) = r q , na ,n b \<l,n a ,n b ) (14) 



,n b \q,n a ,n b ) , (13) can reached in a number of steps less than or equal 

to 2N max (N max + 1), where the amplitudes c rLa . nb are 
nonzero for n a ,n b < A max . 

The algorithm is given by the following sequence of 
operations 



where for n n 



n b 



'0,n a ,n b 
Tl,n a ,n b 



e za cos(6 l /2)co, na .„ i) - 
e~ la cos(0/2)ci,„ O)nb 



ie l/3 sin(6»/2)ci,„ Qi „ b 
- ie~ ll) sia(0/2)co )7lai n fc , 
(15) 

where 9 = fi t and r qt n a ,n b = e^* c qtna , nb for n a -n b ^ n. 
Here we have introduced the phases (a, /?, <p q ) in the tran- 
sition. As discussed below, these phases are needed to 
properly zero out amplitudes during the algorithm. This 
phase control can be implemented by applying short shift 
pulses [26l [33l 00], by an amount Auj q for a duration 
r, providing a controllable phase shift Auj q r. The mi- 
crowave pulse also has a controllable phase S, in terms of 
which a — 4>q — —<p\ — Au! q r/2 and /3 — S + Aw,r/2. 
However, for the states discussed here, such phase control 
is not necessary, and we will assume a = (3 = <f) q = 0. 

By shifting the qubit's frequency into resonance with 
resonator a (uj q — uj a ) or resonator b (u) q — uj b ), the 

interaction terms can be made dominant. Letting A = 
e -iH a t/h^ we let 

-41*)= a q,na.n b \q,n a ,n b ). (16) 

q,n a ,T16 

The operation A performs the mapping 

(lO,n a ,n b = C0 S 6n a C0,n a ,n b ~ * sm #n a c l,n a -l,n b 
ai,n a -l,n b = COs9 na Ci irla ~l,n b ~ 3 Sin 0„ o Co. n a ,n b , 

_ (17) 
where 6 na = y/n^g a t. Similarly, letting B — e t ' H i> t / h j we 
find that 



B\V) = 



b q ,n a ,n b \q,n a ,n b } 



(18) 



performs the mapping 



bo,n a ,n b — COS Qn b C-0,n a ,n b lS\T\-9 nh C\ nanh ^.\ 
bl,n a ,n b -l — COS0n b Cl,n a ,n b -l — iSVO-dn b C0,n a ,n b t 



where 8 nb — y/nt,g b t. Note that these three interactions 
provide all the two-level rotations needed between sets of 
qubit-resonator states. 



U 




(21) 



with 



W„ N a 

u a = II A J R aj> and u hJ = II B ikRb,jk- (22) 



fe=0 



The parameters {t a j, t b ,jk-> Ra,ji an d R b,jk) m this op- 
eration depend on the original amplitudes c v , a<nb and are 
found by solving the inverse evolution: 



3=N h 



|0,0,0). (23) 



As in the state synthesis algorithm of Law and Eberly 
[3H [35], each operation is a two-state rotation (speci- 
fied above). The (inverse) algorithm can then be visual- 
ized as a sequence of transitions on the Fock-state dia- 
gram, each step chosen to zero out a particular element 
\q,n a , n b ) . The solution presented here moves amplitudes 
from right to left column-by-column, from top to bottom 
row-by-row, until finally all of the amplitude is located 
at |0, 0,0). Inverting this sequence produces the desired 
transformation U\Q, 0, 0) — > ^target)- 
In more detail, we have 



rrt — 77t7yt jrt . 

U ~ U a u b .l U b,2 



u, 



b,N b ' 



where 



and 



Ul — R a ,iAiR a 2 A 2 ■ ■ ■ R a N A\ 



(24) 



(25) 



Ki = R hA R hi B n ■ ■ ■ lf js. lf ,s ■ ( 26 ) 

The form of t/J . is chosen to transfer all of the ampli- 
tude of the states in row n b — j to n b = j — 1. For each k, 
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q=0 
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FIG. 1: Schematic set of operations to generate an arbitrary 
state of two coupled resonators. In this Fock state diagram, 
each node corresponds to a specific set of photon numbers 
(n a ,n,b), and contains two levels for the qubit. Interactions 
lead to couplings between these states, indicated by the ar- 
rows between the states. There are three key interactions 
used in this scheme. A transfers quanta between the qubit 
and resonator a by a resonant Jaynes-Cummings-type inter- 
action (solid lines). Similary, B transfers quanta between the 
qubit and resonator b (dashed lines). Finally R (curved ar- 
rows) rotates the qubit for states with n a — rib = n (here with 
n — 0, see text). 



b\ - k transfers the amplitude from |0, k, j) to 1 1, k, j — 1) 



after which R\j^ rotates the state to |0, k,j — 1) 

sequence of operations in * can thus be visualized as 
clearing out row j in the Fock-state diagram, stepping 
from k = N a down to k = 0. This procedure is then 
repeated in U* until all of the amplitude is in row j = 
and thus resonator b is in the ground state (rib — 0) . The 
final sequence of operations in U\ is equivalent to the 
state-synthesis procedure of Law and Eberly [34 for a 
qubit coupled to a single resonator (as modified in 33 ). 
This is illustrated in Fig. [2] for N a = N b = 3. 

An important issue not discussed in [3 6) is the issue 
of relative phases of the various states. In fact, since 
uj a < u q < ujb, all of the states will evolve in time, 
so it is most convenient to specify the algorithm using 
an interaction picture with respect to the Hamiltonian 
H q (0) + H r , as above. However, in the course of the al- 
gorithm it is necessary to zero out the amplitude of a 
particular state by rotating into another state. For ex- 
ample, given a qubit state 



The 



|?A) = cos6»|0) +e" i0 sin( 



1 



(27) 



the amplitude for state |1) can be removed by applying 
the rotation e '^v e -^*/2 r e ^ e -*(4>-f/2W2_ Both 

choices require shifting the relative phases between |0) 
and Fortunately, it is always the case that the two 
states coupled by A, B, and R differ in the qubit excita- 
tion. Thus, one can adjust the relative phases of the two 



FIG. 2: Schematic set of steps in the (inverse of the) state- 
synthesis algorithm. Not all of the Rabi rotations are shown. 



•a 
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1 n a 



FIG. 3: Algorithm to generate the state |*) = |3, 0) + |0, 3} of 
two coupled resonators. The sequence of operations is detailed 
in Table 1. The horizontal, vertical, and curved transitions 
are the interactions A, B, and R (see text). 



relevant states by a short duration shift of the qubit fre- 
quency as described above. These small shifts can there- 
fore be included between each of the various operations 
in U. Specifying these phases is not necessary for the 
following examples. 

The total number of operations in U involves iV a A 
unitaries, (N a + l)N b B unitaries, and N a + (N a + l)N b 
Rabi pulses R. Correspondingly, the maximum total time 
for this algorithm is [54"] 



T max - (N a + N b + N a N b )^r + V^t, 

n f~{ 2 9a\/] 

N„ 

+ (7V a + l)^-^- 



(28) 



This assumes that all of the Rabi and shift pulses are n- 
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TABLE I: Procedure for * = jO, 3, 0) + |0, 0, 3} 



Step 



Parameters 



Quantum State 



Ra,l 




= Tv/2, LU d = L0 


|0,0,0)-» 


1,0,0) 


Ai 


9ata,l 


= tt/2 


|0,0,0>- 


0,1,0) 


Ra,2 


Qtqa,2 


= 7T, Ud = LOi 


|0, 0, 0) + i 


1,1,0) 


A 2 


9ata,2 


= tt/(2v^) 


|0,0,0} + 


0,2,0) 


Ra,3 


Qtqa,3 


= TT, U) d = 0J 2 


|0,0,0>-i 


1,2,0) 


A3 


9a.ta,3 


= vr/(2^3) 


|0,0,0>- 


0,3,0) 


Rt,i 


£Uqb,l 


= TV, UJ d = LV 


-t|l,0,0>- 


0,3,0) 


B x 


9btb,l 


= tv/2 


-|o,o,i>- 


0,3,0) 


Rb,2 


Qtqb,2 


— TV,UJ d = Ui-1 


i|l,0,l) - 


0,3,0) 


B 2 


9btb,2 


= »r/(2v^) 


|0,0,2>- 


0,3,0) 


Rb,3 


fltqb,3 


— TT,iOd = W-2 


-*|1,0,2>- 


0,3,0) 


B z 


9btb,3 


= tt/(2^) 


-|0,0,3>- 


0,3,0) 



See 



pulses. Of course, specific instances of the algorithm can 
have times 
examples. 



have times less than T max , as will be seen in the following 



C. Specific Examples 

The simplest example is the construction of a NOON 
state 

|*noon> = -^|0)®(|iV,0) + |0,A0). (29) 

This state is particularly nice, as the solution for the 
inverse evolution admits the following simplification 



U, 



Rt sn B 



(30) 



'b,j "6,jO J 

The evolution times are chosen to move population from 
|0,0,j) -> |l,0,j-l) -> |0,0, This is done for j = 

N —> 1, after which U\ moves population from |0, N, 0) 
to |0,0,0). Each step is a two-state rotation, and only 
trivial phases appear. We can schematically write this as 
U = (AR) N (BR) N , where all of the rotations have 6 = tt 
except for the first R, which is a 7r/2-pulse. The number 
of steps is 4N, with a time [M] 



NOON 



2N 



1\ TT 



N 



2 J fi ZgaVj 



N 

£ 



(31) 



An explicit list of parameters for the twelve steps needed 
for N — 3 is given in Table [I] 

A more complicated example is given by the maximally 
entangled state 



I*) 



= \n) 8) \N - n). 



(32) 



While this only occupies a sparse region of the Fock-state 
diagram, the algorithm does not have the same simplifi- 
cations or explicit solution as the NOON state. Figure [4] 
shows a numerical solution for the eighteen steps of the 
control sequence for U with = 3. 



D. Superconducting Implementation 

The algorithm described above can be implemented 
by a tunable superconducting qubit, such as the phase 
[IT] or transmon [35] qubit, coupled to two on-chip res- 
onators, as shown in Fig. [5] However, there are a number 
of details that will depend on the specific experimental 
implementation. First, the sequence of operations must 
be "programmed" with a time-dependent control pulse 
for the frequency io q {t) and the microwave drive f2(t). 
Such pulses are shown in Fig. [6] for the NOON state 
sequence with A^ = 3. 

First, when analyzing the fidelity of this control se- 
quence, it is important to note that, for fixed couplings 
g a ,gb, the natural basis to describe the various transi- 
tions is in fact the dressed basis, i.e. those states that 
are true eigenstates of the Hamiltonian for f2 = and Lu q 
equal to some fixed value. While these states are close to 
the product states \q,n a ,nb), this is only an approxima- 
tion. It is convenient to imagine that the couplings g a , gh 
can be turned on and off at the beginning and end of the 
algorithm. This could be achieved by using a tunable 
coupler [43, 44J or by moving the qubit frequency to the 
regime to q > uj a ,^b- 

Second, each of the number-state-dependent Rabi 
transitions will in fact require optimization in both am- 
plitude and frequency for each desired Fock state {n a , rib). 
The amplitude optimization is necessary because of the 
variation of the matrix elements (of a x ) between the 
dressed states, while the frequency optimization is nec- 
essary because the expression for Ud in Eq. (Ill is only 



a leading-order perturbative result for a two-level sys- 
tem. The results of such a partial optimization for the 
control pulse of Fig. [6] is shown in Fig. [7] Here we 
have shown the expectation values (n a ), {rib), and (q) 
(in the uncoupled basis), averaged over a window of 5ns 
to remove high-frequency oscillations (due to fixed cou- 
pling and the dressed basis). Each Rabi pulse increases 
the qubit population, while each swap (A or B) transfers 
the qubit population to the oscillator (thereby increas- 
ing (n a ) or (rih)). This control sequence achieves a state 
equivalent to the A^ = 3 NOON state with a fidelity of 
0.975. Higher fidelity should be possible by using more 
sophisticated Rabi and shift pulses, or optimal control 
techniques. 

Third, the dispersive shifts underlying the number- 
state-dependent transitions for a three-level system are 
in fact quite different than the two- level results. For a 
three-level system with level spacings wqi,wi2, the fre- 
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FIG. 4: Procedure for ^ = X^n=o |0, n, 3 — n). Each block includes the step of the sequence, the parameter values, and the 
relevant transition being performed. In addition, the state populations after the step are shown as a Fock-state diagram, with 
the relevant state highlighted in blue. The parameters were found by numerically solving for the inverse evolution equation. 
The forward sequence reads from top-to-bottom, left-to-right, while the inverse sequence is reversed. 



quency required for |0,n a ,n{,) — > |1, n a ,nb) becomes 

J2 \2 



X n a.,rib) 
J Q1 



n 2 n 2 \' 2 
Wbi + 9a (2n a + 1) + n a 

^01 — U a U! a — LOyi 

+ ^— (2n b + l)+ 9 " X2 n b , (33) 



OJl - UJ b 



UJ h - Ul 12 



with A ~ This dispersive shift depends on the re- 
lationship of both ojqi and W12 with w a and u;;,, with 
somewhat surprising results |29[ 142) . However, many 
of the results described above can be adapted to han- 
dle these complications. By setting the frequencies to 
u a < ^12 < <^oi < it is possible to tune the qubit 
to feel an equal (not opposite) Stark shift from the two 
resonators, so that the resonant condition can be written 
as 



u n = uj q + Acjq + nAtu 



(34) 



with n = n a + rib and Aw < 0. In order to properly 
solve the inverse equations, the only change needed is to 



replace the order of the product in U, 



t 

b,j- 



uL = Rt 



The transitions for row j are now marching along 
columns k from left-to-right, preventing transitions to 
propagate back up from row j — 1 (note that this does not 



affect the NOON state synthesis procedure). Otherwise, 
the algorithm performs quite similarly to that described 
above. 



Resonator A 



Control 
Line 



Qubit 



Resonator B 



FIG. 5: Schematic circuit for generating entanglement be- 
tween two superconducting resonators. Resonator A (blue) 
has a fundamental frequency uj a /2n, while resonator B (red) 
has frequency Ub/2-K. These are each capacitively coupled to 
a tunable qubit (gray) with frequency uu q /2iT, with coupling 
strengths g a and gb- The qubit is controlled by an external 
circuit (black). The theoretical results described in the text 
require uj a < u q < ^b- 
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FIG. 6: Control Sequence for NOON State Synthesis with N 
=3. The upper curve indicates shifts of the qubit frequency 
as a function of time, while the lower indicates the microwave 
pulses applied to the qubit. Here we have set oj q /2n — 7GHz, 
uj a /2Tv = 6.3GHz, u b /2-K = 7.7GHz, and g/2n = 70MHz. The 
various steps of the sequence are labelled, with the first Rabi 
transition a 7r/2-pulse, and the remaining all 7r-pulses. 
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FIG. 7: Expectation values of the qubit excitation (q) and 
photon numbers (n a ) and {rib) as a function of time for the 
NOON state sequence with N = 3, using the control sequence 
and parameters of Fig. [6] These values have been averaged 
over a window of 5ns to remove oscillations due to the dressed- 
state components of the wavefunction (see text). 



III. NOON STATE SYNTHESIS WITH 
DECOHERENCE 

There are two known procedures to generate the en- 
tangled state between two resonators 

|*targct> = ^(|iV,0)+|0,iV)), (36) 

which we shall call Method 1 (for 36J) and Method 2 (for 
[3"2l 133]). The behavior of these two procedures under 
decoherence is the subject of this section. 



We consider the role of dissipation on this process, 
modeled by a Lindbald equation 



dp 
~dt 



■[n, P ] 



3 

£ 



•M /•.;/< - \ L \ L iP~ \p L ] L 3) 



(37) 

where Ai = 1/T 9 , A2 = A3 = l/T r , L\ = er_, L2 = a, 
and L3 = b. This equation is perturbatively solved in the 
Appendix, where the final state of the two resonators has 
the approximate form 

P = Pa,a\N, 0){N, 0| + p a ,b\N, 0)(0, N) 

+ Pbta \0, N){N, 0| + Pbtb \0, N)(0, N\, (38) 

and explicit expressions for the populations p a , a ,Pb,b & n d 
the coherences p a ,b, Pb,a are derived. From these, we cal- 
culate the NOON state fidelity by 



T= (* 



target 



|p|^targct> = \ (Paa + Pab + Pba + Pbb) ■ 



(39) 

Method 1 uses the state synthesis algorithm discussed 
above, for which the preparation sequence can be writ- 
ten as U — {BR 01 ) N {AR 01 ) N . For simplicity, we in- 
clude only resonant interactions for each time-step, as 
described in [29], for both our numerical and analytical 
calculations. As shown in the Appendix, the fidelity can 
be approximated by 



Mi 



exp 

















P.-" 




1] 













x exp 



1 N 

E A *» 



1 2n + N - 1 



T 



T 



,(40) 



where At = ir/Q and At n — tt /(2gy / n). This is shown 
as a function of N in Fig. [8] for typical physical param- 
eters, along with a numerical simulation of the Lindblad 
equation using a quantum trajectories method. Each 
data point was calculated using 1024 trajectories with 
the same physical parameters. The figure demonstrates 
that the perturbative result derived in the Appendix is 
quite good. 

While the dominant loss of fidelity in this method is 
due to dissipation in the qubit, there is still a significant 
contribution related to dissipation int he resonator. The 
quadratic dependence on N of the latter is due to the 
fact that the Fock state \n) decays with a rate n/T r . To 
a good approximation, the overall process has a fidelity 
of the form 



Mi 



a -(7/32)T Ml /T -(l/2)NT Ml /T r 



(41) 



where Tm x = 2NAt + 2 ^ n At n is the total duration of 
the sequence for Method 1. 

The recent NOON state experiment [55] used a differ- 
ent procedure (Method 2), with two three level systems 
initially in the entangled state 



1 



V2 



(|oi) + |io». 



(42) 



and subject to the sequence of pulses of the form 
(AB) [ABR$ Rfz ] JV ~ 1 . Here A, B, and R 12 are detailed 
in the Appendix, where the fidelity for this method is 
shown to be approximately given by 



Method 1 



KI-i 



exp 



x exp 



x exp 



11 



TV — 



N-l 



11 / T, 



n=l 
1 At Q 1 

2^T ~ 2 



2 
2n 



1 



T 



At 



N 



2TV- 1 



(43) 



with At = n/Ct, At n = Tr/{2V2ng) for n = 1 -> JV - 1, 
Aio = 7r/(4g), and At^v = ft/(2>VNg) (see the Appendix 
for more details). This fidelity is shown as a function of 
TV in Fig. [9] using the same parameters as Fig. [8j Here 
the fidelity has the approximate form 



M 2 



-{ll/8)T M JT q -{l/2)NT M JT r . 



(44) 



where Tm2 = NAt + J2 n is the total duration of 
the sequence for Method 2. Note that this method has a 
stronger dependence on the qubit coherence, due to the 
fact that the second-excited state decays at a rate 2/T g . 
The total time of Method 2, however, is approximately 
half that of Method 1, as the (AR)*- 1 and {BR) nJi 
operations are performed in parallel. 

The two procedures, for TV = 4, are directly compared 
in Fig. 10 as a function of the qubit decay time T q . 
For the same values of At = tt/CI, and equal couplings 
g > fl, Method 1 outperforms Method 2 for all values of 
TV and T q . While this is somewhat surprising, it should 
be pointed out that the second method does not use any 
dispersive interactions, so the Rabi pulses can be driven 
faster to decrease the overall preparation time even fur- 
ther. 



IV. CONCLUSION 

In this paper we have explored the synthesis of arbi- 
trary entangled states between two superconducting res- 
onators. Elsewhere, one of us has proposed quantum 
logic operations for such resonators as qudits [29) . Such 
methods, using the larger Hilbert space afforded by har- 
monic oscillator modes, may prove to be an important al- 
ternative to qubit operations. The fabrication of on-chip 
stripline or coplanar waveguide resonators is often much 
easier than qubits based on Josephson junctions, and can 
be expected to yield coherence times of 10/zs, much higher 
than the qubits. Recent experiments |45j have shown 
that three-dimensional cavity resonances can have even 
higher coherence times T r rs 50/iS. Understanding how 
to use the larger resources afforded by these devices re- 
mains a primary goal. We conclude with a discussion of 
important topics for future study of entangled resonators. 

It is interesting to compare and contrast the results of 
the state-synthesis algorithm of [3S] (Method 1) with the 
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FIG. 8: State fidelity T Ml as a function of TV for the NOON 
State Sequence Method 1. The upper curve is for a qubit 
dissipation time T q = 1000ns, while the lower curve is for 
T q — 500ns. Other parameters are resonator decoherence 
T r = 10/is, 0,/2-k = 20MHz, and g/2iv = 100MHz. 



Method 2 
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FIG. 9: State fidelity Jm 2 as a function of TV for the NOON 
State Sequence Method 2. The upper curve is for a qubit 
dissipation time T q = 1000ns, while the lower curve is for 
T q — 500ns. Other parameters are resonator decoherence 
T r = 10/is, Q/2ir = 20MHz, and g/2n = 100MHz. 



NOON state preparation procedure of [3H (Method 
2). Both involve a set of Rabi and swap pulses, and 
exhibit similar results for the fidelity. Method 1, how- 
ever, requires the use of the number-state-dependent 
Rabi transitions in the dispersive regime. These selective 
transitions must have Rabi frequencies less than the sepa- 
ration in frequencies to neighboring transitions, and thus 
high-fidelity operations will typically be slower than non- 
selective Rabi transitions used in method 2. While we 
have not attempted a full optimization of either method, 
it is clear than a faster state synthesis algorithm is desir- 
able, especially in the presence of decoherence. 

Entangled resonators can be used to test higher- 
dimensional Bell inequalities (35] . While the Bell inequal- 
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FIG. 10: State fidelity as a function of qubit dissipation time 
T q for the NOON State Methods with N = 4. The upper 
curve is for Method 1, while the lower curve is for Method 2. 
Other parameters are resonator decoherence T r = 10000ns, 
Q./2-k = 20MHz, and g/2n = 100MHz. 



ity has been tested using two superconducting qubits 
[3TJ27], and the Mermin inequality [35] with three super- 
conducting qubits [40, 49 , there remains much to be ex- 
plored about nonlocality for coupled qudits [50] . Higher- 
dimensional inequalities allow for both more sensitive and 
more robust demonstrations of quantum mechanics, and 
have recently been studied experimentally with optical 
photons [51] for d up to 12. Aside from the entangled 
states themselves, these tests would require qudit logic 
gates [29 and measurements. While qudit measurements 
can be accomplished using coupled qubits in either the 
resonant [33] or quasi-dispersive [39] regimes, these qubit- 
based schemes encode the resonator state in a sequence 
of two-outcome results. It would be advantageous to de- 
velop d-outcome measurements (with d > 2) to more 
efficiently read out the resonator states. 

Entangled resonators may also be useful for quan- 
tum communication. As resonators have the capacity 
to store greater entanglement than qubits, these would 
fit in nicely with the quantum routing protocols of |52j . 
Finally, it may be possible to extend the quantum op- 
tics analogy with circuit-QED further, and use these 
entangled states as interferometric probes to achieve 
Heisenberg-limited measurements [31[ 153) of microwave 
fields. These and other probing questions will undoubt- 
edly lead to new and improved methods for quantum 
control, measurement, and sensing with superconducting 
devices. 
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Appendix: Decoherence Calculations 

Each step of the two NOON state synthesis procedures 
involves a two-state oscillation, where each state is sub- 
ject to decoherence, taken to be dissipation of the qubit 
and oscillator with decay times T q and T r , respectively 
(dephasing can be included similarly). This is modeled 
by a reduced master equation for the density matrix for 
the two states |1) and |2) 

p = p 11 |l)(l|+ Pl2 |l)(2|+p 21 |2)(l|+ /922 |2)(2|. (A.l) 

The master equation takes the form 

dpn 

dt 
dp\2 

dt 
dp2i 



dt 

dp22 

dt 



n 

l 2 {p12 


- P2l) 


- Aipn - 


- Ai 2 p 22 , 


n 

l 2 ipu 


- P22) 


4(Ar4 


A 2 )pi 2 , 


.n 

l^\P22 


- pn) 


4(Ar4 


A 2 )p 2 i, 


.n 


- pn) 


— A 2 p 22 , 





(A.2) 



where £1 is the rotation rate between the two states 
and the various decoherence rates Ai,A 2 ,A i2 are sum- 
marized in Table 1. This equation can be solved ap- 
proximately using standard first-order perturbation the- 
ory (for Ai, A 2 , Ai 2 <C £2) to find 

p n (t) = a + (t) Pll (0)+a^{t)p 22 (0) 

+ib(t) [pi 2 (0)- P21 (0)], 
Pl2 (t) = c+(t)p 12 (0) + c_(t)p 21 (0) 

+ib(t) [p u (0)-p 22 (0)], 
P21 (t) = c+(t)p 21 (0) + c-(t)p 12 (0) 

+ib(t) [p 22 (0)-Ph(0)], 
p 22 (t) = a + (t)p 22 (0) + a-(t)p u (0) 

+ib(t)[p 21 (0)-p 12 (0)], (A.3) 

with 



a±(t) 



-(Ai+A 2 -Ai 2 )*/ 2 



± l e -(2A 1+ 2A 2+ A 12 ) t /4 cos(m); (A4) 

b(t) = 1 e-( 2Al + 2A2 + Al2 ) t / 4 sin(Ot), (A.5) 

C ± (t) = 1 6-^+^/2 

± i e -(2A 1+ 2A 2+ A 12 ) t /4 cos{iU y (A. 6 ) 



Coherences with other states oscillate at frequency fl/2 
and decay with a rate (Ai + A 2 )/4. 
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TABLE II: Parameters for Perturbative Treatment of Deco- 
herence 



after the swap: 



T = p$|00)(00| + p ?|00)<02| +p^|02)(00| +^ ; |02)(02| 

(A.11) 



S 2 ) 



(2) | 



.(2) I 



Operation State 1 State 2 



Ai 



A 2 



Roi |0,n) |l,n> 

i? 12 |l,n> |2,n> 

Ai |0,n+l) \l,n) 



n/T r 1/T, + n/T r 

1/T, + n/T r 2/T, + n/T r 
(n + \)/T T 1/T, + n/T r 



|l,n + l) |2,n> 1/T, + (n + 1)/T r 2/T q +n/T r 



l/2tyith 
2/T, 







1. Method 1 

This method uses a qubit coupled to the two res- 
onators, subject to a sequence of Stark-shifted Rabi tran- 
sitions R i and swapping A,B of excitations from the 
qubit to resonator a and b of the following form 

[BR Q1 (tt)] n [ARo^f- 1 AR 01 (n/2). (A.7) 

Each of the Rabi transitions is a Stark-shifted transition 
resonant for those Fock states with a given n — n a ± rib- 
For simplicity, we assume that all rotations have an equal 
Rabi frequency fi, and that the qubit-resonator couplings 
are also equal g a = 9b = 9- 

We first look at the (AR) N sequence, using the nota- 
tion \q, n a ) for states of the qubit q = 0, 1 and oscillator 
state n a . The first rotation is a 7r/2-pulsc, and generates 
the following state of the qubit-oscillator system. 



p = ^ ) |00)(00|+^|00)(10|+p^|10)(00|+^|10)(10| 

(A.8) 

with = = = pf^ = 1/2. The perturbative 
corrections are in fact negligible for this rotation, and we 
have ignored any overall phases of the coherence terms. 
The excitation is then swapped in time At\ = w/(2g) 
from the qubit to the oscillator, leading to 



(0)| 



p = ^ 1 ) |00)(00|+^|00)(01|+^|01)(00|+piV|01)(01| 

(A.9) 

with 
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(A.10) 



Note that our perturbative approach includes only the 
decaying matrix elements, and not any population of the 
lower excited states. 

The subsequent rotation is a Stark-shifted 7r-pulsc, 
such that |00> |00) and |01) -> |11). The effect of de- 
coherence on this transition is somewhat different, lead- 
ing to the following approximation for the density matrix 
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2 ^T, T r 



PiV- (A.12) 



Here Ai = 7r/f2 is the time for a 7r-pulse, the swap time 
At 2 = Ati/y/2 has decreased, due to the matrix elements 
of the Jaynes-Cummings interaction, and the resonator 
coherence time has also decreased. 

Continuing the sequence of 7r-pulses we find that after 
n steps the state is 



p = p^|00)(00|+p o ^|00)(0n|+p^|0n)(00|+p^|0n)(0n| 

(A.13) 

with the recursion relation 
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(A.14) 



and where Ai„ = Ati/y/n. 

Thus, after the sequence of (Ai?)^ pulses, the final 
qubit-resonator state is 

p = p( r) |00)(00|+pW|00)(07V| 

+P^ ) |0A)(00|+ P ^ ) |0A)(0/Y| (A.15) 
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(A.16) 



We now turn to the (BR) N sequence. Extending the 
notation to \q,n a ,rib), the initial state is 

p = /9 W|000)(000|+pW|000)(0TV0| 

+/4o ) |0M))(000| + p[i ) \0N0){0N0\. (A.17) 

The 7r-pulse is now Stark-shifted so that the only transi- 
tion is 1 000) — > 1 100), or more generally |00n b ) — > |10n 6 ). 
However, the excitations of resonator a continue to de- 
cay, so that the coherences get factors of e - NT / 2T r w hile 
the |iV)(iV| population gets the factor e ~ NT / T -r^ where 
T = A At + J2 n A *« is thc total duration 01 the (BR) N 
sequence. Otherwise the calculation proceeds as above, 
with a final state 

p = p O o\00N)(00N\+ p Q1 \00N)(0N0\ 

+/9io|0M))(00iV| + pn|0A0)(0A0|. (A.18) 



with density matrix elements 



Poo 



1 



Poi 



Pio 



Pn 



x cxp 



exp 



x cxp 



cxp 



x exp 



exp 



-3At/4T, 



N 



N 



cxp 



n=l 



2n- 1 



A V ' T* A ' 



1 N 



T 



■>-!)!- 

n=l x H 



4 T 9 
2n- 1 

4 T 9 
2n- 1 



-N 



(JV) 

Poo > 



2 T r 



(iV) 

Poi ' 



2 T r 



(iV) 

Pio > 



(A.19) 



Removing the qubit state, we have the desired form: 

p = p a , a \N,0)(N,0\+ p atb \N,0)(0,N) 

+p b , a \0, N) (A, 0| + p b , b \0, A) (|0, A|. (A.20) 
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(A.21) 



The fidelity for this state is 

F = (*targot|p|*targot) = ^ {Paa + p a b + Pba + Pbb) ■ 

(A.22) 

Taylor expanding the various terms and regrouping, the 
fidelity can be approximated by 
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(A.23) 



Substituting T = A At + J2 n ^« yields the result stated 
in the text. 



2. Method 2 

In this method, two three-level systems (qutrits) are 
initially prepared in a two-qubit Bell state 

-L(|01> + |10». (A.24) 

Then, each qutrit is rotated from |1) — > |2), followed 
by qubit-resonator swaps (A^ and B^) of the form 
1 2, n) — > |1, n + 1). This sequence can be written as 



A w B w \ a ^b^r^r[ 



N-l 



(A.25) 



12 



where R12 indicates a rotation between the qutrit states 
|1) and |2), and the final swaps (A^ and B^) are of 
the form \l,n) — > |0,n + 1). Note that the swaps and 
rotations can be performed in parallel. For simplicity, we 
assume that all rotations have an equal Rabi frequency Q, 
and that the qubit-qubit and qubit-resonator couplings 
are all equal to g. 

We first observe that the two-qutrit Bell state can be 
formed by a 7r-pulse 1 00) — > 1 10) followed by a partial 
("square-root of") swap 1 10) —> 1 10) + |01) (again ig- 
noring any phases in the various operations). Including 
dissipation, the initial two-qutrit density matrix is ap- 
proximately 



P 
2 

where 



(|01)(01) + |01)(10| + |10)(01| + |10)(10|) (A.26) 
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-3At/4T, 



-Ato/T, 



(A.27) 



with At — tt/Q and At = 7r/(4g). 

Extending the notation to \q a , qi,, n a , n&), the first R12 
pulse yields the state 



p = p$ |0200) (0200| +p$ |0200) (2000| 
+ / o^ ) |2000)<020( 
with density matrix elements 



+pfj 1 2000) (0200 1 +p ( $ 1 2000) (2000 1 (A.28) 
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(A.29) 



This is then swapped to the resonators, so that the total 
state becomes 

p = P w 10101) (01011 +/$ 10101) (10101 

-f-p^l 1010) (0101 1 +^ 1 1 ) l 1010 >( 101 °l ( A - 30 ) 

with 



where At\ = n/(2^/2g). The factor of \^2 in Ati is due 
to the qutrit matrix element (l|cr_|2) |29) . 

After n steps of this sequence, the state becomes 

p = p { ^\0l0n)(010n\+ p { ^\010n)(l0n0\ 

+p ( $ |10ra0) (010n| + p$ | lOnO) (10n0| (A.32) 

where the density matrix elements satisfy the recursion 
relation 



(n+l) 

Poo 



J.n+1) 
P01 



Pio +1) 



J.n+1) 
Pll 



1 e -At/2T, 



2 

x exp 



1 

2 e 



-2At/T q 



x e 



At 1l+ i 



T, 



e -3At/2T qe -nAt/T r 



x exp 



At„ H 



T 



e -3At/2T qe -nAt/T r 



x exp 
1 

2*" 
x exp 



At n+ i 



2n + 1 



2n + 1 
T 

r 



2n + 1 



-nAt/T r 



(n) 

Poo » 



(n) 
Pol 



(n) 
PlO > 



At/2T q 



At n 



i e -2At/T g \ e -nAi/T r 



+ 1 



2n + 1 



P11 » 
(A.33) 



and Ai„ = n/(2\/2ng). This sequence repeats until n = 
N-l. 

The final step swaps the remaining qubit excitations 
into the resonators, yielding the state 

p = p ( Q o ) \000N}(000N\+ p { o ^\000N)(00N0\ 

+p[q } I00A0) (000A| + p[^ |00iV0) (00A0| 

(A.34) 

with 



P01 



pfo } 



Pll 



P ( 1 -At/2T q 

2 V 2 



x exp 



E e -3At/2T q exp 
P_ e ~3At/2T q exp 
P ( 1 -At/2T, _ 

2 V 2 



-e 
2 



-2At/T q 



-2At/T q 





Aii 






x exp 


"~2~ 







Ati 






-2- 






Ati 






-"2- 




4). 



(A.31) 



(AT) 

Poo 



(AT) 

P01 



(AT) 
PlO 



Pll 



exp 
exp 
exp 
exp 



At 



N 



2 

At N 
2 

At N 
2 

At N 



1 2N-1 



T, 



1 2N-1 



1 2N -I 



T 



Tr 



1 2N-1 



T„ 



Tr 



{N-l) 

Poo 1 



(N-l) 

P01 > 



{N-l) 
PlO . 



Pll ' 



(A.35) 



and Aijy = n/feyNg). This leaves the qutrits in their 
ground states, and tracking the various decoherence fac- 
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tors wc find 



The fidelity for this state is 



P ( 1 p -M/2T q , 1 p -2M/T q 



Pab 



Pba 



2 \2 

x exp 

x exp 
P 

2 ex P 

x exp 
P 

2 ex P 
x exp 



N-l 



1 N 

-;2> 



3 — 28 n ,,\ 2n — 1 



1 w 



3 - 25 n . N 2n-l 
T T 

n -*- r 



1 N 



n=l 



3 — 2<5„ : 7v 2n — 1 



2 \2 
x exp 

x exp 



iV-1 



■ i(iV-l)(JV-2)^ 



n=l v 



3 — 2S n ^N 2n — 1 



(A.36) 



-T 7 = (*targot|/3|*targot> = ^ (Paa + Pab + Pba + Pbb) ■ 

(A.37) 

Taylor expanding the various terms and regrouping, the 
fidelity can be approximated by 



T 



exp 



x exp 



x exp 



11/ T„ 



N-l 



n=l 

1 At 



2 

2n- 1 



1 



,1 2AT-1 

2 r„ -2 A ^It; + ^- 



(A.38) 
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